Quantum limited spin diffusion in the unitary Fermi gas 
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We compute spin transport in the unitary Fermi gas using the strong-couphng Luttinger-Ward 
theory. In the quantum degenerate regime the spin difTusivity attains a minimum value of Ds — 
1.3h/m approaching the quantum limit of diffusion for a particle of mass m. Conversely, the spin 
drag rate reaches a maximum value of Fad — 1.2 ksTp /h in terms of the Fermi temperature Tf- 
The frequency-dependent spin conductivity (Ts(i^) exhibits a broad Drude peak, with spectral weight 
transferred to a universal high-frequency tail (Ts(w — > oo) = C /3Tv(muj)^^^ proportional to the Tan 
contact density C. For the spin susceptibility Xs{T) we find no downturn in the normal phase above 
Tc ~ 0.16 Tf. 
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The excitation and decay of spin currents plays an im- 
portant role in many areas of condensed matter physics, 
including the development of electronic devices based on 
spin transport. While the Coulomb interaction does not 
affect electrical currents in a uniform system [1] it trans- 
fers momentum between spin-up and down particles and 
thereby dampens the spin current. Understanding the 
mechanism of spin drag and spin diffusion quantitatively 
is important for an effective control of spin currents, how- 
ever in solids this is often complicated by the presence of 
impurities and lattice effects. Ultracold atomic Fermi 
gases provide an extremely clean experimental realiza- 
tion to study the effect of the two-particle interaction 
alone If the interactions are short-ranged and the 
scattering length is much larger than the particle spac- 
ing the results are universal and apply to a wide range of 
models, including even dilute nuclear matter. The spin 
diffusivity Dg measures how quickly a spin current levels 
out a gradient in the spin density. In a strongly interact- 
ing Fermi gas Dg decreases as the temperature is lowered 
into the quantum degenerate regime and reaches a min- 
imum near the Fermi temperature Tp, before increasing 
again at even lower temperatures in the superfluid phase. 
The minimum value of Ds in the strong-coupling region 
can be understood qualitatively as a consequency of the 
uncertainty principle: the mean-free path in a gas cannot 
become shorter than the mean particle spacing in the ab- 
sence of localization |3l| , which translates into a quantum 
bound Dg }^ Ti/m for particles of mass to. For a strongly 
interacting Fermi gas of trapped ^Li atoms a spin diffu- 
sivity Ds > 6.3(3) ?i/m has recently been measured [1]. 
Very low spin diffusion is found also in graphene [5J , while 
spin Coulomb drag in GaAs quantum wells yields a value 
of Ds > 500 h/m [g. 

The determination of Ds near its minimum in the 
strongly interacting regime, and more generally the ques- 
tion whether quantum mechanics imposes universal lower 
bounds on the transport coefficients, is a key challenge in 
many-body physics. Recent progress comes from the Anti 



de Sitter and Conformal Field Theory correspondence 
which maps a strongly coupled field theory to an equiv- 
alent weakly coupled gravitational theory, where calcu- 
lations are feasible. It gives a lower quantum bound for 
the internal friction of mass flow, expressed as the ratio 
of shear viscosity to entropy ry/s > h/AnkB, in certain 
relativistic field theories [7[. Quantum limited friction, 
or perfect fluidity [1], has been found to be almost sat- 
isfied in very different physical situations ranging from 
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quark-gluon plasmas to ultracold atomic gases 
remains an open question whether a similar bound ex- 
ists for spin diffusion in non-relativistic systems [l^ . In 
this Letter we present a strong-coupling calculation of 
the spin diffusivity Ds in the homogeneous unitary Fermi 
gas, which is the only known example of a non-relativistic 
interacting scale-invariant fluid. In Fig. [1] our result is 



shown by the solid red line, and we find a minimum value 
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FIG. 1: (Color online) Spin diffusivity Ds (in units of Ti/m) 
versus reduced temperature T/Tf (solid red line). The ex- 
perimental data 0] (blue squares) for the trapped gas are 
rescaled down by a factor of 4.7 to compensate for the ef- 
fect of the trapping potential. The dashed black line is 
the result for a classical gas of weakly interacting particles, 
Ds = 1.1 {T/TFf^^h/m. 
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of Dg ~ \.?>h/m at a temperature of about T = 0.5 Tp. 
To our knowledge this is the lowest value achieved to date 
for a strong two-particle interaction. Recent experimen- 
tal data for the trapped unitary gas 3 are shown as the 
blue squares. The data are rescaled down by a factor of 
4.7 to compensate for the effect of the trapping potential, 
and we obtain a remarkable agreement for all tempera- 
tures in the normal phase above the superfluid transition 
temperature ~ 0.16 Tp- 

In the high-temperature limit where Ds 3> h/m is 
much larger than the quantum limit our calculation 
agrees with the predictions of Boltzmann kinetic theory 
13-31 (dashed black line). In the strongly interact- 



ing region near Tc, however, the fermions cease to be 
well-defined quasi-particles 0, [l3| and Boltzmann the- 
ory is not applicable. Therefore, we employ the strong 
coupling Luttinger-Ward theory to compute spin trans- 
port. The Luttinger-Ward formalism 18, 13] accurately 
describes both the normal and the superfluid phase of 
the BEC-BCS crossover problem [g^: the values for 
Tc/Tf = 0.16(1) and the Bertsch parameter ^ = 0.36(1) 
agree within error bounds with recent precision experi- 
mental [111 and diagrammatic Monte Carlo [22| results. 
By construction it fulfills the conservation laws including 



scale invariance f9| and the Tan relations exactly |11|. 
The Luttinger-Ward theory has recently been extended 
to compute transport coefficients in linear response us- 
ing the Kubo formula: this gives access to the frequency- 
dependent shear viscosity of the unitary Fermi gas, which 
was found to satisfy the exact viscosity sum rule 0, . 
We now extend this work to the case of spin transport 
in order to explain the recent experiment by Sommer 
et al. 4], and we proceed as follows: first we compute 
the frequency-dependent spin conductivity (7s(a;) of the 
unitary Fermi gas. The dc value CTs = crs(w = 0) de- 
termines the spin drag rate T^d = n/mas at density n, 
which is the rate of momentum transfer between atoms 
of opposite spin. We then compute the spin suscepti- 
bility Xs — <9(n-[- — ni)/d{^i^ ~ /ij.) which characterizes 
the magnetic properties of the system 13, 13] ■ Finally, 
we determine the spin diffusivity shown in Fig. [T] by the 
Einstein relation Dg ~ o's/Xs- 

The strongly interacting two-component Fermi gas is 
described by the grand canonical Hamiltonian 
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V 



fc,fc' ,q 

(1) 



where Sk — /2m (h = 1) is the free particle disper- 
sion and Ha the chemical potential for the cr = t, | com- 
ponents. The s-wave contact interaction 170 acts only 
between different fermion species at low temperatures. 
The bare interaction 170 is singular in the ultraviolet 
and needs to be regularized; the renormalized coupling 
g = Anh^a/m determines the s-wave scattering length a. 



The transport coefficients are obtained from the mi- 
croscopic model via the retarded number-current/spin- 
current correlation function 



Xjn/js(g,w) 



dt / d^xe'^'^'~i-'''> 



[j^±jl]{x,t),[j^±jl]{O,0) ) . (2) 



The spin selective current operators in Fourier represen- 
tation are given by j^iq) = T,ki^/'^)(i^q/2,a^k+q/2.a- 
The correlation function determines the conductivity 
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ImXj„/js(q,w) 



(3) 



which measures the relaxation of a global number/spin 
current at frequency u. The total response integrated 
over all frequencies is proportional to the particle density 
by the number/spin /-sum rule (25I. [26l| 



duj n 

cr„/s(cj) = — 



(4) 



For a momentum-conserving interaction ([TJ the particle 
current cannot decay and (Tn(w) = ■Kn5{Lj) / m. In con- 
trast, scattering transfers momentum between f and \. 
particles so that the spin current relaxes and (Js{^) has 
a nontrivial structure. 

We compute the current correlation function ([2]) us- 
ing field theoretical methods and Feynman diagrams 
in the Matsubara formalism p5| . The current oper- 
ator = ± j| implies a current response vertex 

Jaa' = Jaa' + + ^aa' the Fevuman diagrams 
which splits into three contributions 0, [3| (c, a' are 
the spin indices of incoming and outgoing fermion lines) . 
The first term is the bare number (spin) current vertex 
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iJ^a,gip) = p.4a') with the £ 
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partial wave component of the momentum p and Pauli 
matrices . The other two terms are current vertex 
corrections which are required to fulfill the conservation 
laws exactly. The Maki-Thompson (MT) contribution 
describes direct scattering between quasi-particles while 
the Aslamazov-Larkin (AL) term captures the induced 
current of fermion pairs, or molecules (for details see 
Ref. i). For a mass current both f and i fermions 
move in the same direction and induce a current of pairs, 
leading to a sizeable AL term. In contrast, for a spin 
current f and X atoms move in opposite directions 
and no pair current is induced. Hence, the Aslamazov- 
Larkin correction to the spin current vanishes exactly 
in the spin balanced case, J^a's — 0, which constitutes 
an important simplification. We solve the self-consistent 
equation for the fully dressed current vertex Ja-a' by it- 
eration and obtain the current correlation function ([2|) 
via the Kubo formula Since the correlation func- 
tion Xjn/js('7 = 0,iWm) is evaluated at discrete imaginary 
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FIG. 2: (Color online) Spin conductivity (Ts(a;) (in units of 
Tin/mEp) versus frequency (red circles) at T = 0.5 Tf. The 
Drude model (solid black line) has the same total spectral 
weight as (Ta(cj) given by the spin /-sum rule. Part of the 
spectral weight is transferred from lower frequencies into a 
universal high-frequency tail (dotted blue line) (7s(aj oo) = 
C/37r(maj)^''^ with Tan contact density C = 0.086 [gj. 



FIG. 3: (Color online) Spin drag rate Fsd (in units of Ef/K) 
versus reduced temperature T /Tf (solid red line). The ex- 
perimental data (blue squares) for a trapped gas are 
rescaled up by a factor of 5.3 to compensate for the ef- 
fect of the trapping potential. The dashed black line is 
the result for a classical gas of weakly interacting particles, 
Fsd = 0.9(r/I»-i''2£F/ft. 



Matsubara frequencies we must perform an analytic 
continuation in order to obtain the physically relevant 
correlation function Xjn/js(w) for real frequencies w. We 
use Pade approximants and find that the continuation 
is robust at low temperatures if we vary the number of 
Matsubara frequencies. In addition, it yields the correct 
high-frequency tail (see below). Specifically, we oversam- 
ple the Matsubara data twice with a spline fit and use 
the first five Matsubara frequencies in order to extract 
the spin drag rate Fsd- We validate our strong coupling 
calculation by confirming that crs(w) indeed fulfills the 
spin /-sum rule Q within 1%. 

Spin conductivity. The resulting spin conductivity 
CTs(w) is shown in Fig.[2]for reduced temperature T /Tp = 
0.5 where it has the lowest dc value cts = 0.8 n/m (red cir- 
cles). In a Drude model the conductivity would assume a 
form af™'^^(w) = (n/m)Fsd/(w2 + F^^) (solid black line) 
with total spectral weight given by the sum rule. The 
spin drag rate Fgd is a parameter which we determine 
from the dc limit CTs — n/mTsd of our full numerical so- 
lution. We find that the true as{uj) deviates from the 
Drude model for lo }^ Ep: spectral weight is transferred 
from the region uj 8 Ep to higher frequencies where it 

(dotted blue 



forms a power-law tail as{uj — > oo) 
line in Fig.©. 

The high-frequency response generally depends on the 
non-universal short-distance behavior of the interatomic 
potential. However, for a broad Feshbach resonance as in 
^Li this potential has a range much shorter than the 
particle spacing, fcF|?'e| ^ 1, and becomes effectively a 
contact interaction. In this case the correlation functions 
exhibit universal power-law tails in the high-frequency 
range ma.x{Ep,kBT) < i:j < h/{mrl) [27] which de- 



pend only on the Tan contact density C 28|. In the 
high-frequency limit the exact transport equations can 
be solved analytically in a manner analogous to the vis- 
cosity response [0| , and we obtain the universal spin con- 
ductivity tail 



as{ijJ —5- oo) 



C 



37r(mw)-^/2 



(5) 



in agreement with the result from the operator product 
expansion (29| . Similar tails appear in other transport 
properties such as the viscosity 0, 2^ 2^, 3^ . The value 
for the Tan contact density C = 0.0863 k% at T/Tp = 0.5 
extracted from the tail of (Ts(w) agrees better than 1% 
with the value C = 0.0860 from the tail of the mo- 
mentum distribution ~ Ck~'^ [^. A similar behavior 
of (Js{ijj) is observed for all temperatures T >Tc- 

We now turn to the dc limit and plot the spin drag 
rate Fgd = n/mcTs in Fig. [3] (solid red line). The spin 
drag has a maximum value of Vgd ~ 1.2 Ep/h in the 
quantum degenerate regime around T /Tp = 0.5 and de- 
creases both for lower and higher temperatures. In the 
high-temperature limit of a classical gas the Luttinger- 
Ward transport equations can be solved analytically to 
leading order in the fugacity and we obtain Fgd = 
(32V2/97r3/2)(T/TF)-i/2^^/;i = 0.9 (T/Tp)-^/^ Ep/h 
for T ^ Tp in agreement with Boltzmann kinetic theory 
0, [13] • The fact that the numerical solution at large tem- 
peratures agrees with the analytical result for T ^ Tp 
is a nontrivial validation of our analytical continuation 
procedure. 

The measured spin drag rate in a trapped unitary 
Fermi gas Q (blue squares in Fig. has the same quali- 
tative behavior as our numerical data, with a broad maxi- 
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mum between T/Tp = 0.4 . . . 0.8. Note that the absolute 
spin drag rate cannot be directly compared to our calcu- 
lation for the uniform system: the solution of the trans- 
port equation depends on the trap geometry and the ve- 
locity profiles of f and | particles in the trap [i , 1^ 31 1 . 
For a quadratic velocity profile in a harmonic trap the 
spin drag rate T'^^^^ = Tsd/a is rescaled by a constant 
factor a = 2^/^ in the high-teniperature limit (see sup- 
plementary information of Ref . [J] ) . In the experiment a 
factor of a = 5.6(4) is found, and we obtain the best fit at 
high temperatures for a = 5.3. In the quantum degener- 
ate regime T ^ Tp the assumption of a uniform quadratic 
velocity profile breaks down: in the center of the trap a 
large spin drag leads to slow spin motion, while the spins 
in the weakly interacting wings move rapidly. The ve- 
locity profile thus becomes non-uniform and a acquires 
a temperature dependence. In Fig. [3] the calculation for 
the uniform system and the trap averaged data rescaled 
by a constant factor differ for T ^ Tp. We thus find 
that the scaling factor starts to deviate from the high- 
temperature estimate a = 5.3 around Tp. 

Spin susceptibility. We shall compute and discuss the 
spin susceptibility Xs in order to find Dg = ag/xs- Both 
the spin susceptibility Xs = d{n-f — ni)/d{fi-f — fii) and 
the normalized compressibility Xn = n^K = d(ni^ + 
ni)/d{fj,-^ + fj^l) are obtained from the number/spin cor- 
relation function 



Xn/s 



dtd^a 



[nf±ni]{x,t),[n^±ni]{O,0) 



(6) 



The spin-selective particle number operator in Fourier 
representation is given by n„{q) = J2k ^k~q/2.a^k+q/2,a- 
In the Luttinger-Ward formulation the bare number 
(spin) density vertex has the form N^„,^ = r°^, {N^„,^ = 
T^^,). For the dressed number density vertex both 
MT and AL vertex corrections contribute, while the 
AL term again vanishes for Ng in the spin balanced 
case. In order to obtain the static susceptibility Xn/s 
we calculate the susceptibility Xn/s{i^m) for zero Mat- 
subara frequency, i.e., Xn/s = Xn/s(*Wm = 0). Note 
that an analytical continuation is not needed here. The 
static limit of the related current correlation function 
Xjn/js = Xjn/js(«'^m = 0) = n/m IS fixcd by the exact /- 
sum rule Q , and our numerical computation fulfills this 
sum rule within 1%. We therefore expect our results for 
the static susceptibilities Xn/s to be of the same accuracy. 

The susceptibilities of the free Fermi gas are Xn = Xs = 
n/T for T Tp (Curie- Weiss) and x° = Xs = Xo = 
3n/2Ep for T ^ in the Fermi liquid phase (dashed 
black line in Fig. |4]). In the unitary Fermi gas the at- 
tractive interaction leads to a compressiblity Xn twice as 
large as Xn in the quantum degenerate regime near T^ 
(open magenta circles), in very good agreement with the 
experimental data [4] (full cyan circles) and with a non- 
self-consistent diagrammatic approach [3^ . Conversely, 
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FIG. 4: (Color online) Compressibility Xn (circles) and spin 
susceptibility Xs (squares) versus reduced temperature T/Tp. 
The experimental data [3| (full symbols) is compared to 
our Luttinger-Ward calculation (open symbols). The dashed 
black line is the susceptibility of the free Fermi gas. 



we find that the spin susceptibility Xs remains below Xs 
and exhibits a maximum of about Xs — O-dn/Ep at 
T/Tp = 0.3 (open red squares). The spin susceptibility 
is expected to vanish as exp(— 2A/r) deep in the super- 
fluid phase with gap A, where an infinitesimal magnetic 
field gradient cannot break pairs. The proposed pseudo- 
gap scenario 17 , S^l predicts a pronounced drop of Xs at 
a pair- breaking scale T* > T^. Our data, which fully in- 



clude the attractive branch, remain nearly constant down 
to T/Tp ~ 0.2. This indicates that the scales T* and T^ 
are very close in the unitary Fermi gas. The measured 
Xs (full blue squares in Fig. 3]) also shows no downturn 
and can be described in a Fermi liquid picture despite the 
large value for T^/Tp ~ 0.16, in agreement with a recent 
quantum Monte Carlo and experimental study 33] . Note 
that a finite condensate fraction can lead to significantly 



lower values for Xs [34 1 



In the experiment by Sommer et al. Q Xs is determined 
from a combination of the local spin density gradinet 
and the trap-averaged center of mass motion. Hence, the 
measured Xs deviates at low temperatures from the cal- 
culation of the uniform system (cf. Fig. 2]) for similar rea- 
sons as discussed above for Fgd. Remarkably, we find that 
the differences in Tgd and Xs cancel and lead to a very 
good agreement in the spin diffusivity Dg = n/mFgdXs 
shown in Fig. [1] above. 

In conclusion, our strong coupling calculation of spin 
transport explains the behavior of the spin diffusivity Dg 
seen in experiment and we find that the diffusivity 
of the unitary Fermi gas reaches the quantum limit. We 
predict a univeral high-frequency tail of the spin conduc- 
tivity (Ts(w) which should be accessible experimentally 
using Bragg spectroscopy for the dynamic structure fac- 



tor 35| . It would be desirable to have local measurements 
of transport properties in a way similar to local precision 
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measurements of the thermodynamic properties 21| and 
the momentum distribution [361] . 
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